clc;
[x,y,z] = meshgrid(0:0.1:2);
V = sqrt(x.^2 + y.^2 + z.^2);
slice(x,y,z,V,[1 2],[1 2],[0 1]);

% ---------------------------------------------------------------------------
[x0,y0]=meshgrid(0:0.1:2);  % 生成一个从0到2，步长为0.1的二维网格数据
z0=ones(size(x0));          % 生成一个与x0同样大小的全1矩阵
h=surf(x0,y0,z0);          % 绘制一个表面图，初始为z=1的平面
rotate(h,[1,0,0],45);       % 沿x轴正向旋转45度

x1=get(h,'XData');         % 提取该平面数据的x坐标
y1=get(h,'YData');         % 提取该平面数据的y坐标
z1=get(h,'ZData');         % 提取该平面数据的z坐标

slice(x,y,z,V,x1,y1,z1),hold on,slice(x,y,z,V,2,2,0)
[x,y,z]=meshgrid(0:0.1:2);
V=sqrt(x.^x+y.^((x+y)/2)+z.^((x+y+z)/3));
vol_visual4d(x,y,z,V)